Single-cell transcriptome profiling reveals cell type-specific variation and development in HLA expression of human skin.

Skin, the largest organ of body, is a highly immunogenic tissue with a diverse collection of immune cells. Highly polymorphic human leukocyte antigen (HLA) molecules have a central role in coordinating immune responses as recognition molecules. Nevertheless, HLA gene expression patterns among diverse cell types within a specific organ, like the skin, have yet to be thoroughly investigated, with stromal cells attracting much less attention than immune cells. To illustrate HLA expression profiles across different cell types in the skin, we performed single-cell RNA sequencing (scRNA-seq) analyses on skin datasets, covering adult and fetal skin, and hair follicles as the skin appendages. We revealed the variation in HLA expression between different skin populations by examining normal adult skin datasets. Moreover, we evaluated the potential immunogenicity of multiple skin populations based on the expression of classical HLA class I genes, which were well represented in all cell types. Furthermore, we generated scRNA-seq data of developing skin from fetuses of 15 post conception weeks (PCW), 17 PCW, and 22 PCW, delineating the dynamic expression of HLA genes with cell type-dependent variation among various cell types during development. Notably, the pseudotime trajectory analysis unraveled the significant variance in HLA genes during the evolution of vascular endothelial cells. Moreover, we uncovered the immune-privileged properties of hair follicles at single-cell resolution. Our study presents a comprehensive single-cell transcriptomic landscape of HLA genes in the skin, which provides new insights into variation in HLA molecules and offers a clue for allogeneic skin transplantation.

Biomarkers of Innate Immunity and Immunological Susceptibility to Viral Infection in Patients with Alcoholic Cirrhosis.

BACKGROUND: The harmful effect of alcohol on the immune system may be due to both a direct action of the alcohol or its metabolites on immune cells as an indirect action modifying the different mechanisms of intercellular interaction. The interplay between stimulatory (aKIR) and inhibitory (iKIR) natural killer (NK) cell receptors and their corresponding human leukocyte antigen (HLA) ligands influences the outcome of virus infection. The aim was to analyze the influence of the KIR/HLA pair genetic profile in male alcoholic cirrhosis (AC) patients with and without viral infections to find susceptibility biomarkers that can help establish the risks and prevent viral infections. METHODS: A total of 281 male AC patients were analyzed. The sociodemographic characteristics, viral hepatitis C (HCV), hepatitis B (HBV), and cytomegalovirus (CMV) infections were analyzed. Genomic DNA was extracted, and genetic the KIR/HLA profiles were investigated. A total of 6 KIR genes and their corresponding ligands (HLA-C) were analyzed. Patients were grouped into two groups: with and without associated viral infection. RESULTS: A statistically significant increase in the combination of KIR2DL2+/C1C1 was observed in male AC patients with viral infection compared to those without viral infection (45.9% vs. 24.5%, p = 0.021). The analysis of KIR2DL3+/C1+ showed a high frequency comparing healthy controls and male AC patients without virus infection (85% vs. 76.4%; p = 0.026). The analysis of KIR2DL3+/C2C2 frequency showed a statistically significant increase comparing male AC patients without viral infection and healthy controls (23.6% vs. 15%; p = 0.026). CONCLUSIONS: The genetic KIR2DL2+/C2C2 profiles may play a significant role in determining the vulnerability of male AC patients to viral infections, providing valuable insights for future research and potential therapeutic interventions.

The gut microbiome molecular mimicry piece in the multiple sclerosis puzzle.

The etiological complexity of multiple sclerosis, an immune-mediated, neurodegenerative disease with multifactorial etiology is still elusive because of an incomplete understanding of the complex synergy between contributing factors such as genetic susceptibility and aberrant immune response. Recently, the disease phenotypes have also been shown to be associated with dysbiosis of the gut microbiome, a dynamic reservoir of billions of microbes, their proteins and metabolites capable of mimicring the autoantigens. Microbial factors could potentially trigger the neuroinflammation and symptoms of MS. In this perspective article, we discussed how microbial molecules resulting from a leaky gut might mimic a host's autoantigen, potentially contributing to the disease disequilibrium. It further highlights the importance of targeting the gut microbiome for alternate therapeutic options for the treatment of MS.

Single-Nucleotide Polymorphisms Within Non-HLA Regions Are Associated With Engraftment Effectiveness for Patients With Unrelated Cord Blood Transplantation.

Clinically, stem cells with matched human leukocyte antigens (HLAs) must be selected for allogeneic transplantation to avoid graft rejection. However, adverse reactions still occur after cord blood transplantation (CBT). It was inferred that the HLA system is not the only regulatory factor that may influence CBT outcomes. Therefore, we plan to investigate whether the single-nucleotide polymorphisms (SNPs) located in non-HLA genes are associated with the effectiveness of CBT. In this study, the samples of 65 donors from CBT cases were collected for testing. DNA sequencing was focused on the SNPs of non-HLA genes, cytotoxic T-lymphocyte-associated protein 4 (CTLA4), CD28, tumor necrosis factor ligand superfamily 4 (TNFSF4), and programmed cell death protein 1 (PDCD1), which were selected in regard to the literatures published in 2017 and 2018, which indicated that they were related to stem cell transplantation. Then, in combination with the detailed follow-up transplantation tracking database, these SNPs were analyzed with the risk of mortality, relapse, cytomegalovirus (CMV) infection, and graft-versus-host disease (GVHD). We found that there were 2 SNPs of CTLA4, 1 SNP of TNFSF4, and 2 SNPs of PDCD1 associated with the effectiveness of unrelated CBT. These statistically significant SNPs and haplotypes would be used in clinical to choose the best donor for the patient receiving CBT. Moreover, the polygenic risk scores (PRSs) with these SNPs could be used to predict the risk of CBT adverse reactions with an area under the receiver operating characteristic curve (AUC) of 0.7692. Furthermore, these SNPs were associated with several immune-related diseases or cancer susceptibility, which implied that SNPs play an important role in immune regulation.

Construction of an HLA Classifier for Early Diagnosis, Prognosis, and Recognition of Immunosuppression in Sepsis by Multiple Transcriptome Datasets.

Background: Sepsis is a clinical syndrome, due to a dysregulated inflammatory response to infection. Accumulating evidence shows that human leukocyte antigen (HLA) genes play a key role in the immune responses to sepsis. Nevertheless, the effects of HLA genes in sepsis have still not been comprehensively understood. Methods: A systematical search was performed in the Gene Expression Omnibus (GEO) and ArrayExpress databases from inception to 10 September 2021. Random forest (RF) and modified Lasso penalized regression were conducted to identify hub genes in multi-transcriptome data, thus we constructed a prediction model, namely the HLA classifier. ArrayExpress databases, as external validation, were utilized to evaluate its diagnostic, prognostic, and predictive performance. Immune cell infiltration score was calculated via CIBERSORTx tools and single-sample gene set enrichment analysis (ssGSEA). Gene set variation analysis (GSVA) and ssGSEA were conducted to determine the pathways that are significantly enriched in different subgroups. Next, we systematically correlated the HLA classifier with immunological characteristics from multiple perspectives, such as immune-related cell infiltration, pivotal molecular pathways, and cytokine expression. Finally, quantitative real-time polymerase chain reaction (qRT-PCR) was conducted to validate the expression level of HLA genes in clinical samples. Results: A total of nine datasets comprising 1,251 patients were included. Based on RF and modified Lasso penalized regression in multi-transcriptome datasets, five HLA genes (B2M, HLA-DQA1, HLA-DPA1, TAP1, and TAP2) were identified as hub genes, which were used to construct an HLA classifier. In the discovery cohort, the HLA classifier exhibited superior diagnostic value (AUC = 0.997) and performed better in predicting mortality (AUC = 0.716) than clinical characteristics or endotypes. Encouragingly, similar results were observed in the ArrayExpress databases. In the E-MTAB-7581 dataset, the use of hydrocortisone in the HLA high-risk subgroup (OR: 2.84, 95% CI 1.07-7.57, p = 0.037) was associated with increased risk of mortality, but not in the HLA low-risk subgroup. Additionally, immune infiltration analysis by CIBERSORTx and ssGSEA revealed that B cells, activated dendritic cells, NK cells, T helper cells, and infiltrating lymphocytes (ILs) were significantly richer in HLA low-risk phenotypes, while Tregs and myeloid-derived suppressor cells (MDSCs) were more abundant in HLA high-risk phenotypes. The HLA classifier was significantly negatively correlated with B cells, activated dendritic cells, NK cells, T helper cells, and ILs, yet was significantly positively correlated with Tregs and MDSCs. Subsequently, molecular pathways analysis uncovered that cytokine-cytokine receptor (CCR) interaction, human leukocyte antigen (HLA), and antigen-presenting cell (APC) co-stimulation were significantly enriched in HLA low-risk endotypes, which was significantly negatively correlated with the HLA classifier in multi-transcriptome data. Finally, the expression levels of several cytokines (IL-10, IFNG, TNF) were significantly different between the HLA subgroups, and the ratio of IL-10/TNF was significantly positively correlated with HLA score in multi-transcriptome data. Results of qRT-PCR validated the higher expression level of B2M as well as lower expression level of HLA-DQA1, HLA-DPA1, TAP1, and TAP2 in sepsis samples compared to control sample. Conclusion: Based on five HLA genes, a diagnostic and prognostic model, namely the HLA classifier, was established, which is closely correlated with responses to hydrocortisone and immunosuppression status and might facilitate personalized counseling for specific therapy.

How HLA diversity is apportioned: influence of selection and relevance to transplantation.

In his 1972 paper 'The apportionment of human diversity', Lewontin showed that, when averaged over loci, genetic diversity is predominantly attributable to differences among individuals within populations. However, selection can alter the apportionment of diversity of specific genes or genomic regions. We examine genetic diversity at the human leucocyte antigen (HLA) loci, located within the major histocompatibility complex (MHC) region. HLA genes code for proteins that are critical to adaptive immunity and are well-documented targets of balancing selection. The single-nucleotide polymorphisms (SNPs) within HLA genes show strong signatures of balancing selection on large timescales and are broadly shared among populations, displaying low FST values. However, when we analyse haplotypes defined by these SNPs (which define 'HLA alleles'), we find marked differences in frequencies between geographic regions. These differences are not reflected in the FST values because of the extreme polymorphism at HLA loci, illustrating challenges in interpreting FST. Differences in the frequency of HLA alleles among geographic regions are relevant to bone-marrow transplantation, which requires genetic identity at HLA loci between patient and donor. We discuss the case of Brazil's bone marrow registry, where a deficit of enrolled volunteers with African ancestry reduces the chance of finding donors for individuals with an MHC region of African ancestry. This article is part of the theme issue 'Celebrating 50 years since Lewontin's apportionment of human diversity'.

Epidemiology and genetics of granulomatosis with polyangiitis.

Granulomatosis with polyangiitis (GPA) previously known as Wegener's granulomatosis (WG) is a rare rheumatic disease affecting subjects of all ages. Prevalence and incidence of this systemic disease greatly varies across different ethnic groups. GPA is the commonest form of ANCA-associated vasculitis (AAV) with PR3 positivity among 85-95% of the cases. Scientific investigations of GPA is warranted because its severity, clinical heterogeneity, fast disease manifestation and end-organ damage. The etiology of GPA is still unknown. Major role of HLA and non-HLA genes with immune functions were identified, however, very limited replication was observed in different ethnic populations. In the present review, we have discussed the updates on the global epidemiology and contribution of HLA and major non-HLA genes/loci in GPA. We have also highlighted the cross disease association of GPA associated genes that may help in better disease management and predictive medicine. We proposed that high-resolution HLA typing and development of genetic risk model would help in early disease diagnosis and understanding the prognosis.

Relevance of HLA gene polymorphisms in Romanian patients with chronic renal insufficiency undergoing renal transplantation.

BACKGROUND: Chronic renal insufficiency (CRI) is a global public health problem with a high incidence in the Romanian population. In this study, we aimed to investigate genomic HLA polymorphisms in Romanian patients with CRI waiting for kidney transplantation. To determine the existence of a potential strong link between certain HLA polymorphisms and CRI, we also looked at HLA specificity combinations within the same locus or even different loci, referring to randomly inherited allelic combinations rather than potential haplotypes. METHODS: A total of 2199 patients with CRI on the kidney transplantation waiting list were included. A total of 2786 healthy individuals were included as controls. Both patients and controls were assessed for both HLA I and class II genes. HLA genes were typed using the low-resolution method polymerase chain reaction sequence-specific primer. RESULTS: Certain class I and class II HLA allele groups, genotypes and haplotypes were significantly more frequent in patients with CRI than in the control individuals (eg B* 40 (p ≤ .001, pc ≤ .001), C* 12 (p ≤ .001, pc ≤ .001), DRB1*14 (p = .0022, pc = .04), C*12,- (p < .001, pc < .001), A*01-C*15 (p = .0003, pc = .03) and A*02-C*12 (p = .0005, pc = .0486)). CONCLUSIONS: HLA gene polymorphisms could be clinically relevant CRI-associated genetic profiles in Romanian patients with CRI.

TNF-α and α-synuclein fibrils differently regulate human astrocyte immune reactivity and impair mitochondrial respiration.

Here, we examine the cellular changes triggered by tumor necrosis factor alpha (TNF-α) and different alpha-synuclein (αSYN) species in astrocytes derived from induced pluripotent stem cells. Human astrocytes treated with TNF-α display a strong reactive pro-inflammatory phenotype with upregulation of pro-inflammatory gene networks, activation of the nuclear factor κB (NF-κB) pathway, and release of pro-inflammatory cytokines, whereas those treated with high-molecular-weight αSYN fibrils acquire a reactive antigen (cross)-presenting phenotype with upregulation of major histocompatibility complex (MHC) genes and increased human leukocyte antigen (HLA) molecules at the cell surface. Surprisingly, the cell surface location of MHC proteins is abrogated by larger F110 fibrillar polymorphs, despite the upregulation of MHC genes. Interestingly, TNF-α and αSYN fibrils compete to drive the astrocyte immune reactive response. The astrocyte immune responses are accompanied by an impaired mitochondrial respiration, which is exacerbated in Parkinson's disease (PD) astrocytes. Our data provide evidence for astrocytic involvement in PD pathogenesis and reveal their complex immune reactive responses to exogenous stressors.

Pharmacogenetic and Association Studies on the Influence of HLA Alleles and Rivastigmine on the Iranian Patients with Late-Onset Alzheimer's Disease.

Alzheimer's disease (AD) is a progressive neurodegenerative disorder affecting cognitive function. A number of allelic genes from HLA complex have shown variable associations with AD in different populations. In this study, we investigated the association of DQB1*06:00/x, DRB1*04:00/x, DRB1*15:00/x, and B*07:00/x genotypes with AD and their relevance to the efficacy of rivastigmine treatment in the Iranian population. Our findings suggest that DQB1*06:00/x genotype offers strong protection against AD (P = 0.0074), while B*07:00/x genotype imposes a significant susceptibility for sporadic Alzheimer's disease (SAD) (P = 0.009). Interestingly, B*07:00/x genotype does not show any apparent associations with familial Alzheimer's disease (FAD). Our studies also suggest a pharmacogenetic relationship between drug treatment and presence of a particular genotype in the Iranian LOAD patient population. The Clinical Dementia Rating analysis showed that LOAD patients carrying DRB1*04:00/x genotype tend to display a downward trend in the disease severity and symptoms after 2-year follow-up with rivastigmine treatment. Moreover, in our total patient population, the carriers of DQB1*06:00/x and B*07:00/x alleles have better and worse responses to rivastigmine respectively. We also measured the clinical relevance of the testing for these genotypes employing prevalence-corrected positive predictive value (PcPPV) formula. The PcPPV of testing for DQB1*06:00/x in the Iranian LOAD patients was 1.17% which means that people carrying this genotype have half of the probability of the absolute risk for developing LOAD, whereas the PcPPV of testing for B*07:00/x was 4.45% for SAD, which can be interpreted as a doubling chance for developing LOAD among the Iranian population carrying this genotype. These results also suggest that DQß1 peptide containing positively charged AAs histidine30 and arginine55 and HLA class I ß chain containing negatively charges aspartic acid114 and glutamic acid45,152 in their binding groove plays important roles in protection against and susceptibility for LOAD respectively.

Identical twins with idiopathic membranous nephropathy.

The pathogenesis of idiopathic membranous nephropathy is associated with autoantibodies, most often against the phospholipase A2 receptor (PLA2R) and with genetic factors, especially those involving human leukocyte antigen (HLA) genes. Idiopathic membranous nephropathy is not a typical inherited Mendelian disorder. Reports of idiopathic membranous nephropathy in twins are rare. Herein, we report on two twin sisters diagnosed with PLA2R-associated idiopathic membranous nephropathy. We identified the HLA-DRB1*0301, HLA-DRB1*1501, and HLA-DQB1*0602 alleles in the twin sisters, which were reported as independent risk alleles for idiopathic membranous nephropathy in the Asian population. This case report provides novel evidence for the role of predisposing HLA alleles in the pathogenesis of idiopathic membranous nephropathy.

Beyond the HLA Genes in Gluten-Related Disorders.

Most common food grains contain gluten proteins and can cause adverse medical conditions generally known as gluten-related disorders. Celiac disease is an immune-mediated enteropathy triggered by gluten in individuals carrying a specific genetic make-up. The presence of the human leukocyte antigens (HLA)-DQ2 and HLA-DQ8 haplotypes together with gluten intake is a necessary, although not sufficient, condition, to develop celiac disease. Fine mapping of the human genome has revealed numerous genetic variants important in the development of this disease. Most of the genetic variants are small nucleotide polymorphisms located within promoters and transcriptional enhancer sequences. Their importance is underlined by an increased risk in DQ2/DQ8 carriers who also have these non-HLA alleles. In addition, several immune-mediated diseases share susceptibility loci with celiac disease, shedding light on the reasons for co-occurrence between these diseases. Finally, most of the genes potentially involved in celiac disease by fine genetic mapping of non-HLA loci were confirmed in gene expression studies. In contrast to celiac disease, very little is known about the genetic make-up of non-celiac wheat sensitivity (NCWS), a clinically defined pathology that shares symptoms and gluten dependence with the celiac disease. We recently identified differentially expressed genes and miRNAs in the intestinal mucosa of these patients. Remarkably, the differentially expressed genes were long non-coding RNAs possibly involved in the regulation of cell functions. Thus, we can speculate that important aspects of these diseases depend on alteration of regulatory genetic circuits. Furthermore, our finding suggests that innate immune response is involved in the pathogenic mechanism of NCWS. This review is intended to convey the idea that in order to fully understand celiac disease and its relationship with other gluten-related disorders, it is worth learning more about non-HLA variants.

Diagnostic Test Criteria for HLA Genotyping to Prevent Drug Hypersensitivity Reactions: A Systematic Review of Actionable HLA Recommendations in CPIC and DPWG Guidelines.

INTRODUCTION: Certain HLA variants are associated with an increased risk of hypersensitivity reactions to specific drugs. Both the Clinical Pharmacogenetics Implementation Consortium (CPIC) and the Dutch Pharmacogenetics Working Group (DPWG) have issued actionable HLA gene - drug interaction guidelines but diagnostic test criteria remain largely unknown. We present an overview of the diagnostic test criteria of the actionable HLA - drug pairs. METHODS: A systematic literature search was conducted in PubMed, Embase, Web of Science and Cochrane Library. Original case-control and cohort studies were selected and sensitivity, specificity, positive predictive value (PPV), negative predictive value (NPV) and number needed to genotype (NNG) were calculated for the actionable HLA-drug pairs. RESULTS: In general, the HLA tests show high specificity and NPV for predicting hypersensitivity reactions. The sensitivity of HLA tests shows a wide range, from 0-33% for HLA-B*1502 testing to predict lamotrigine induced SJS/TEN up to 100% for HLA-B*5701 to predict immunologically confirmed abacavir hypersensitivity syndrome (ABC-HSR). PPV is low for all tests except for HLA-B*5701 and ABC-HSR which is approximately 50%. HLA-B*5701 to predict ABC-HSR shows the lowest NNG followed by HLA-B*5801 for allopurinol induced severe cutaneous adverse drug reactions and HLA-B*1502 for carbamazepine induced SJS/TEN. DISCUSSION: This is the first overview of diagnostic test criteria for actionable HLA-drug pairs. Studies researching HLA genes and hypersensitivity are scarce for some of the HLA-drug pairs in some populations and patient numbers in studies are small. Therefore, more research is necessary to calculate the diagnostic test criteria more accurately.

Psoriatic arthritis: A systematic review of non-HLA genetic studies and important signaling pathways.

Psoriatic arthritis (PsA) is a common, chronic inflammatory disease with complex pathogenesis. In recent years, a number of susceptibility non-human leukocyte antigen (HLA) genes of PsA have been revealed, which also act as important factors in the pathogenesis of PsA as well as HLA genes. By searching the databases National Center for Biotechnology Information, Google and PubMed, 37 articles are included and 50 susceptibility non-HLA genes for PsA are presented, such as IL23A, TNIP1, TYK2, STAT4, IL12B, RUNX3 and TRAF3IP2. In these non-HLA genes, some are common genes shared with other diseases, whereas most of these susceptibility genes are related to the pathogenesis of PsA by activation or inhibition of the signaling pathways. Several signaling pathways possibly implicated in the pathogenesis of PsA are introduced in this paper, including the 2 mainly signaling pathways, IL23/Th17 signaling pathway and NF-κB signaling pathway, and the other involved signaling pathways, such as JAK-STAT signaling pathway and MAPK signaling pathway.

A new set of reagents and related software used for NGS based classical and non-classical HLA typing showing evidence for a greater HLA haplotype diversity.

To evaluate the HLA typing performance of a new Long-Range PCR NGS set of reagents and its dedicated software, a panel of 41 reference homozygous cell lines from the International Histocompatibility Working Group (IHWG) and a panel of 376 volunteer bone marrow donors were analyzed for classical and non-classical HLA class I and class II genes. All results, except HLA-DPB1, were obtained without any ambiguities at the 3rd field level. Based on the high resolution performance of the reagents, a number of new alleles have been described not only for classical but also for non-classical HLA class I genes, leading to a more accurate haplotype definition. Linkage disequilibrium between HLA-A and HLA-G genes has been defined at 4th field level of resolution. Moreover, for the first time, HLA-DQA2 and DQB2 polymorphisms and their linkage disequilibrium with DQB1 were described.

A Computational Platform and Guide for Acceleration of Novel Medicines and Personalized Medicine.

In the era of big data and informatics, computational integration of data across the hierarchical structures of human biology enables discovery of new druggable targets of disease and new mode of action of a drug. We present herein a computational framework and guide of integrating drug targets, gene expression data, transcription factors, and prior knowledge of protein interactions to computationally construct the signaling network (mode of action) of a drug. In a similar manner, a disease network is constructed using its disease targets. And then, drug candidates are computationally prioritized by computationally ranking the closeness between a disease network and a drug's signaling network. Furthermore, we describe the use of the most perturbed HLA genes to assess the safety risk for immune-mediated adverse reactions such as Stevens-Johnson syndrome/toxic epidermal necrolysis.

Human Leukocyte Antigens (HLA) Genes Association in Type 1 Diabetic Nephropathy.

BACKGROUND: Diabetic nephropathy is a common worldwide multifactorial disease where involvement of genetic factors is well etablished. The aim of this study was to investigate the HLA genes implication in the development of type 1 diabetic nephropathy. METHODS: We performed a case- control study where one hundred and fifty subjects were examined. Patients were divided in two groups; with and without type 1 diabetic nephropathy. HLA typing was performed using Polymerase Chain Reaction- Sequence Specific Oligonucleotide (PCR- SSO) method. HLA association to clinical phenotype and HLA haplotype analysis was also investigated. RESULTS: HLA B*51 is increased in patients without type 1 diabetic nephropathy (7.14% vs. 0 %, P <0.05, OR= 0), however no other studied alleles seem to have any effect (all P>0.05). Haplotype analysis also does not reveal any significant association, however, A*02-B*18-DRB1*03-DQA1*05- DQB1*03 haplotype shows a tendency to be associated with the development of diabetic nephropathy (P = 0.05). CONCLUSION: These results suggest a protective effect of HLA B*51 allele from type 1 diabetic nephropathy. However, further studies are required in order to clarify its potential implication as a protective marker.

Influence of Genetic Polymorphism Towards Pulmonary Tuberculosis Susceptibility.

Tuberculosis (TB) is still remains the major threat for human health worldwide. Several case-control, candidate-gene, family studies and genome-wide association studies (GWAS) suggested the association of host genetic factors to TB susceptibility or resistance in various ethnic populations. Moreover, these factors modulate the host immune responses to tuberculosis. Studies have reported genetic markers to predict TB development in human leukocyte antigen (HLA) and non-HLA genes like killer immunoglobulin-like receptor (KIR), toll-like receptors (TLRs), cytokine/chemokines and their receptors, vitamin D receptor (VDR) and SLC11A1 etc. Highly polymorphic HLA loci may influence antigen presentation specificities by modifying peptide binding motifs. The recent meta-analysis studies revealed the association of several HLA alleles in particular class II HLA-DRB1 with TB susceptibility and valuable marker for disease development especially in Asian populations. Case-control studies have found the association of HLA-DR2 in some populations, but not in other populations, this could be due to an ethnic specific association of gene variants. Recently, GWAS conducted in case-control and family based studies in Russia, Chinese Han, Morocco, Uganda and Tanzania revealed the association of genes such as ASAP1, Alkylglycerol monooxygenase (AGMO), Forkhead BoxP1 (FOXP1), C-terminal domain phosphatase 1 (UBLCP1) and intergenic SNP rs932347C/T with TB. Whereas, SNP rs10956514A/G were not associated with TB in western Chinese Han and Tibetan population. In this review, we summarize the recent findings of genetic variants with susceptibility/resistance to TB.